Solvable model for chimera states of coupled oscillators 



OO 
O 

o 

(N 

< 

(N 

Q 

U 
d 



(N 
> 

in 
o 

o 

OO 

o 



Daniel M. Abrams, 1 Renato E. Mirollo, 2 Steven H. Strogatz, 3 'El and Daniel A. Wiley 4 

'Department of Earth, Atmospheric, and Planetary Sciences, 54-621, 
Massachusetts Institute of Technology, Cambridge, MA 02139, USA. 
2 Department of Mathematics, Boston College, Chestnut Hill, MA 024 67, USA. 
'Department of Theoretical and Applied Mechanics, Cornell University, Ithaca, NY 14853, USA. 
^Department of Mathematics, University of Maryland, College Park, MD 20742, USA. 

(Dated: August 26, 2008) 

Networks of identical, symmetrically coupled oscillators can spontaneously split into synchronized 
and desynchronized sub-populations. Such chimera states were discovered in 2002, but are not well 
understood theoretically. Here we obtain the first exact results about the stability, dynamics, and 
bifurcations of chimera states by analyzing a minimal model consisting of two interacting populations 
of oscillators. Along with a completely synchronous state, the system displays stable chimeras, 
breathing chimeras, and saddle-node, Hopf and homoclinic bifurcations of chimeras. 
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Many creatures sleep with only half their brain at a 
time Such unihemispheric sleep was first reported in 
dolphins and other sea mammals, and has now been seen 
in birds and inferred in lizards When brain waves 
are recorded, the awake side of the brain shows desyn- 
chronized electrical activity, corresponding to millions of 
neurons oscillating out of phase, whereas the sleeping side 
is highly synchronized. 

From a physicist's perspective, unihemispheric sleep 
suggests the following (admittedly, extremely idealized) 
problem: What's the simplest system of two oscillator 
populations, loosely analogous to the two hemispheres, 
such that one synchronizes while the other does not? 

Our work in this direction was motivated by a se- 
ries of recent findings in nonlinear dynamics [!, 0, 0, 
0, 0, Hi- In 2002, Kuramoto and Battogtokh reported 
that arrays of nonlocally coupled oscillators could sponta- 
neously split into synchronized and desynchronized sub- 
populations The existence of such "chimera states" 
came as a surprise, given that the oscillators were iden- 
tical and symmetrically coupled. On a one-dimensional 
ring 0, 0] the chimera took the form of synchronized do- 
main next to a desynchronized one. In two dimensions, 
it appeared as a strange new kind of spiral wave [H , with 
phase-locked oscillators in its arms coexisting with phase- 
randomized oscillators in its core — a circumstance made 
possible only by the nonlocality of the coupling. These 
phenomena were unprecedented in studies of pattern for- 
mation 0] and synchronization (Io| in physics, chemistry, 
and biology, and remain poorly understood. 

Previous mathematical studies of chimera states have 

assumed that they are statistically stationary 1,0,0,0, 

0] • What has been lacking is an analysis of their dynam- 
ics, stability, and bifurcations. 

In this Letter we obtain the first such results by con- 
sidering the simplest model that supports chimera states: 
a pair of oscillator populations in which each oscillator 
is coupled equally to all the others in its group, and less 



strongly to those in the other group. For this model 
we solve for the stationary chimeras and delineate where 
they exist in parameter space. An unexpected finding is 
that chimeras need not be stationary. They can breathe. 
Then the phase coherence in the desynchronized popula- 
tion waxes and wanes, while the phase difference between 
the two populations begins to wobble. 

The governing equations for the model are 
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where a = 1,2 and N a is the number of oscillators in 
population a. The oscillators are assumed identical, so 
the frequency lo and phase lag a are the same for all of 
them. The strength of the coupling from oscillators in a' 
onto those in a is given by K aa i . To facilitate comparison 
with earlier work, we suppose that K\\ = K22 = ^ > 0, 
and K\2 = K21 = v > 0, with fi > v. Thus, the coupling 
within a group is stronger than the coupling between 
groups. This corresponds to the assumption 0, 0, 0] 
of a nonlocal coupling that decreases with distance. By 
rescaling time, we may set fx + v = 1. It also proves 
useful to define the parameters A = fx — v and (3 = ir/2 — 
a, because, as we'll show, chimeras exist only if these 
quantities are small enough. 

Simulations of Eq. (TT]) display two types of behavior. 
For many initial conditions, the system approaches the 
synchronized state where all 8's are equal. Otherwise it 
evolves to a chimera state (Fig. [1}. The oscillators in 
group 1 arc in sync; those in group 2 are not. 

Figure [2] illustrates the dynamics of chimera states. 
We plot the phase coherence of the desynchronized pop- 
ulation, as quantified by the order parameter r(t) = 
|^ e ^j(*)) 2 | j w here the angle brackets denote an average 
over all oscillators in population a = 2. In Fig. [2]Ja) 
the order parameter remains constant, except for slight 
fluctuations due to finite-size effects. Thus, this chimera 
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FIG. 1: Snapshot of a chimera state, obtained by numerical 
integration of CTJ with (3 = 0.1, A = 0.2, and Ni = iV 2 = 
1024. The two populations are shown side by side. 



is stable and statistically stationary. However, if we in- 
crease fi (the coupling within a population) relative to v 
(the coupling between populations), the stationary state 
can lose stability. Now the order parameter pulsates, and 
the chimera starts to breathe (Fig. [2Kb)). The breathing 
cycle lengthens as we increase the disparity A = n — v 
between the couplings (Fig. HJc)). At a critical disparity, 
the breathing period becomes infinite. Beyond that, the 
chimera disappears and the synchronized state becomes 
a global attractor. 

To explain these results, we analyze Eq. {l} in the con- 
tinuum limit where N a — > oo for a = 1,2. Then Eq. ([1]) 
gives rise to the continuity equations 
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where f 7 (9,t) is the probability density of oscillators in 
population er, and v a (9,t) is their velocity, given by 



v a (0,t) =uj+ 



2 



\(6'-6-a)f a '{6',t)d0'. (3) 



(Note that we dropped the superscripts on 9 to ease the 
notation. Thus, 9 means 9 a and 9' means 9 a .) If we 
define a complex order parameter 
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then v a {9) simplifies to 

v a (0,t) = w + lm[e 
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FIG. 2: Order parameter r versus time. In all three panels, 
^1 = ^2 = 128 and f3 = 0.1. (a) A = 0.2: stable chimera; 
(b) A — 0.28: breathing chimera; (c) A — 0.35: long-period 
breather. Numerical integration began from an initial condi- 
tion close to the chimera state, and plots shown begin after 
allowing a transient time of 2000 units. 



where the * denotes complex conjugate 
Following Ott and Antonsen 



11] 



we now consider a 
special class of density functions f a that have the form 
of a Poisson kernel. The remarkable fact that Ott and 
Antonsen discovered is that such kernels satisfy the gov- 
erning equations exactly, if a certain low-dimensional sys- 
tem of ODEs is satisfied. In other words, for this family of 
densities, the dynamics reduce from infinite dimensional 
to finite (and low) dimensional. (Numerical evidence sug- 
gests that all attractors lie in this family, but proving this 
remains an open problem.) 
Specifically, let 
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What's special here is that we use the same function a a (t) 
in all the Fourier harmonics, except that a a is raised to 
the n th power in the n th harmonic. Inserting this f a into 
the governing equations, one finds that this is an exact 
solution, as long as 



0. 
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Instead of infinitely many amplitude equations, we have 
just one. (It's the same equation for all n.) 

To close the system, we express the complex order pa- 
rameter z a in terms of a a . Inserting the Poisson ker- 
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nel ((6|) into Eq. ([4]), and performing the integrals, yields 



r (t)= ^ a*At), 
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by orthogonality. Thus the amplitude equations become 
1 



= di + iu)a\ + —a 2 (Knal + K^a^) e 
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and similarly for &2, with l's and 2's interchanged. 

Rewrite the amplitude equations using polar coordi- 
nates p and (j>, defined by a a = Pa&^ % ^" , a = 1,2. (The 
negative sign is included in the definition of <f> so that the 
Poisson kernel f a converges to 5(9 — <p a ), not 5(6 + (p a ), 
as p — > 1 from below. Thus <p a can be interpreted as the 
"center" of the density f , and p a measures how sharply 
peaked it is.) Then Eq. © becomes 
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= -pi4>! + PiUJ 

1 + p 2 

r— [/xpi sin a + ^2 sin (<p\ — (p 2 + a)] (10) 

and similarly for p2 and 4>2, with l's and 2's interchanged. 

Now consider the case of chimera states, for which one 
population is in sync while the other is not. Taking a = 1 
to be the synchronized population, we set pi = 1, corre- 
sponding to a (5-function for that group. Note that p\ = 1 
satisfies the governing equations for all time, since p\ = 
when pi = 1. Hence the condition pi = 1 defines an in- 
variant manifold, on which the dynamics reduce to 

0i = u> — p sin a — vr sin(^) + a) 
1 - r 2 



[/zr cos a + ^ cos("0 — a)] 
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2r 



[/xr sin a + i/ sin(o! — V*)] (11) 



where we've defined r = p2 and tp = <p\ — <p2- The (r, ^>) 
dynamics decouple, yielding a 2-£> system given by 
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[//r cos a-\- v cos(-0 — a)] 
[/xr sin a. — v sin(0? — a)] — p sin a 
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This system has a trivial fixed point r = 1 , tp = 
(perfectly synchronized state) which always exists. The 
non-trivial fixed points correspond to stationary chimera 
states, in which the local order parameters p\(t) = 1 and 
P2(t) = r(i) remain constant, as docs the phase difference 
ip(t) = <pi(t) — 02 (t) , despite the fact that the individual 
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FIG. 3: Phase portraits for Eq. (|12p . regarding r and tp as po- 
lar coordinates. Parameters as in Figs. E|a),(b), respectively, 
(a) Stable chimera (solid diamond), (b) Breathing chimera, 
shown as a stable limit cycle (thick curve) about unstable 
chimera (open dot). In both panels: open diamond, saddle 
chimera; thin solid line, unstable manifold; dashed line, stable 
manifold; solid dot, stable synchronized state. 



microscopic oscillators in population a = 2 continue to 
move in a desynchronized fashion. 

Figure [3] plots typical phase portraits for (fT2j) . Fig- 
ure E^a) shows a stable chimera state coexisting with 
the stable synchronized state; the basin boundary be- 
tween them is defined by the stable manifold of a saddle 
chimera. As we increase the disparity A between the cou- 
plings p and v, the chimera becomes less stable and even- 
tually undergoes a supercritical Hopf bifurcation, creat- 
ing a stable limit cycle (Fig.[3Jb)), the counterpart of the 
breathing chimera of Fig. 2(b). 

Additional phase plane analysis (not shown) reveals 
two other bifurcations. With further increases in A, the 
limit cycle expands and approaches the saddle. Mean- 
while, its period lengthens, which accounts for the be- 
havior seen earlier in Fig. [UJc). At sufficiently large A 
the cycle touches the saddle point and destroys itself in 
a homoclinic bifurcation. On the other hand, if A is de- 
creased from its value in Fig.[3ja), the stable and saddle 
chimeras in Fig.[3fa) approach each other, and eventually 
coalesce and annihilate in a saddle-node bifurcation. 

Figure 0] summarizes the bifurcations and stability re- 
gions for the system. In the rest of the paper we outline 
the analysis leading to these results. 

To calculate the fixed points for Eq. (fT2"|) , we set r = 
and r ^ 1 (since group a = 2 is desynchronized) and 
obtain pr cos a + ^cos(a — tp) = 0. Substituting p = 
(1 + A)/2, v = (1 - A)/2 and (3 = tt/2 - a into r = and 
solving for A yields 



A 



sin(/3 + ip) + r sin (3 
sin(/3 + tp) — r sin/3 



(13) 



at a fixed point. Then imposing tp = and using the 
expression for A above, we find 



sin(2/3 + tp) 



1.(2/3 + 2 sin V>' 



(14) 
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FIG. 4: Stability diagram for chimera states. Bifurcation 
curves: saddle-node (dotted) and supercritical Hopf (solid), 
both found analytically; homoclinic (dashed), found numeri- 
cally. All three curves intersect at a Takens-Bogdanov point 

(P,A) = (cos- 1 V VT3/8+i/2,^^^l) = (0.2239,0.3372). 



Equations (fT3|) and (fT4|) together parametrize all the 
fixed points. They define a surface, or equivalcntly, a 
two-parameter family. Sweeping /3 and ip yields the cor- 
responding r and A values. The resulting surface has two 
sheets that join along a fold; its projection onto the (/3, A) 
plane defines the curve of saddle-node bifurcations. 

To calculate the saddle-node curve, we linearize (fT2")) 
about a fixed point, and set the determinant of the Ja- 
cobian to zero. This implies 

. sin(2/3 + ip) [sin(/3 - 2ip) + 2 sin(/3 + 2ip)] n 

sm 3 H ; = 0, 

' sin(2/3- ip) + 2sin?/> 

(15) 

where we've used Eqs. (fT3|) and (jT4j) to simplify the de- 
terminant. Solving (fT5|) for VK/^) yields two roots, but 
one of them implies r > 1 and hence is spurious; the 
correct root is 



-2/3-2/3 2 + 2/3 3 



11/3 4 



3271/3 6 
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12/3 5 - ""'^ + 0(/3 7 ). 

(16) 

This is then substituted into (fT4"|) to yield r(/3), which 
in turn yields A(/3), via (|13[) . In this way we obtain the 
saddle-node curve 
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-0(/3 6 ), (17) 



which matches the numerical curve shown in Fig. 21 

To find the Hopf curve, we set the trace of the Jacobian 
to zero, which gives ip = — \ sin -1 (2 sin 2/3). Repeating 
the procedure above leads to an exact parametric equa- 



tion for the Hopf curve. Its leading order behavior is 



Ag.{fi) = 2-\/3+(4V3 - 6) (3 2 + (-^= - lO^j /3 4 +0(/3 6 ). 

(18) 

Future work should investigate whether breathing 
chimeras exist for the one- and two-dimensional ar- 
rays of oscillators studied previously 0, Q, d, @, 0]. 
Are the stability diagrams for such systems similar to 
Fig. 4? Do chimeras also exist if the oscillators are 
non-identical [l(], EH [ll] or arranged in complex net- 
works |13|? It would also be worth looking for experi- 
mental examples of chimera states. Candidate systems 
include arrays of lasers 14 1 and chemical [l5| or electro- 
chemical [16j oscillators. 
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